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Abstract 

We demonstrate that the two-dimensonal electron system in a strong per- 
pendicular magnetic field has stable states which break rotational but not 
translational symmetry. The Laughlin fluid becomes unstable to these states 
in quantum wells whose thickness exceeds a critical value which depends on 
the electron density. The order parameter at 1/3 reduced density resembles 
that of a nematic liquid crystal, in that a residual two-fold rotation axis is 
present in the low symmetry phase. At filling factors 1/5 and 1/7, there 
are states with four- and six-fold well. We discuss the experimental 

detection of these phases. 
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The construction of the Laughlin fluid and its identification with the fractional quantum 
Hall effect (FQHE) marked the discovery of a qualitatively new many-body phase. [I] Until 
that discovery, it was generally assumed that the only phases present in the two-dimensional 
electron system were the usual liquid phase and the Wigner solid, both of which are present 
in the classical one-component Coulomb plasma phase diagram. A hexatic phase is also a 
possibility in this system as well as in the logarithmic potential case perhaps more relevant 
to the FQHE. || This theoretical background, and the experimental discovery of phases 
which seem to have properties unlike both the Laughlin liquid and the Wigner solid, @] 
("Hall insulator") led us to an investigation of the possibility of broken rotational symmetry 
(BRS) in the two-dimensional electron system in strong field. We found that several liquid 
crystal-like phases occur, though these appear to resemble more closely nematic than hexatic 
phases. 

We begin with the case where the electron density n satisfies n = §2^72, wnere ^ is the 
magnetic length: £~ 2 = eB/hc and B is the external field, taken to be in the negative z- 
direction. In this paper we work in the limit of large field. This is the density for the 1/3 
quantum Hall state. Consider the following wavefunction for the disc-shaped system: 

N 

= U.[{zi - Zj){zi - Zj - ai){ Zi - z 3 + a£)) exp(- ]T \zi\ 2 /i£ 2 ). (1) 

i<j i 

Here Zj = Xi + iyi and i,j are particle indices, a is a complex number. This wavefunction 
is antisymmetric in the particle indices for all a, lies entirely in the lowest Landau level, 
and reduces to the Laughlin wavefunction at a = 0. It also shares with the Laughlin 
wavefunction the characteristic of having uniform density (far from the edges of the disc). 
However, the two-particle correlation function is a different matter: 

„ _ N(N-l)Ol*l 2 fd 2 z t )\^ 

In this equation f = {x\ — x%, y\ — y^), and N is the total number of electrons . Translational 
symmetry is not broken for the small values of a with which we are concerned, as we shall 
demonstrate below. Thus g a is a function only of the difference variable r. In contrast to the 



Laughlin liquid, however, it does depend on the direction of r. Let us take a to be real. Then 
the equivalent classical plasma interaction corresponding to ^> a is a logarithmic interaction 
between rodlike charged objects lying along the x-axis. Accordingly, represents a BRS 
state. H In general (Rea, Ima) is a director, not a vector, order parameter, since "$? a = ^f- a . 
Thus g a {f) does not have full rotation symmetry for a ^ but does satisfy g a (x,y) = 
9a(-x,y) = g a (x, -y) = g a (-x, -y). 

To give a physical picture of this state, we display a typical configuration in a Monte 
Carlo simulation governed by the probability distribution \^ a \ 2 in Fig. [l], with a = 3.2. One 
notes immediately the stripes along the director, reminiscent of a nematic state. This value 
of a is unphysically large and is chosen for purposes of illustration only. 

Under what conditions is such a BRS state stable? Or, given an interelecron potential 
V{fi — fj), when do we have U a = (^ a \V\^ a ) < U = (^ \V\^o) ? It is easy to see, 
expanding the polynomial in Eqn. [I], and taking the limit \r\ — > in Eqn. |2] that g a {r) ~ r 2 
rather than g a (r) ~ r 6 . Thus for very short-range V, the Laughlin state is always favored, 
as is well known. |J [0 

In order to compare the a = with the a ^ states we performed Monte Carlo 
simulations of the equivalent classical plasmas with 200 particles in the disk geometry, and 
compared the energies of candidate ground states. (We believe that improvement of the 
wavefunction by, for example, quantum Monte Carlo techniques |§ will not affect energy 
differences very much.) In order to minimize the finite size effect we compute the energies 
of only 25 particles closest to the center of the disk. We have checked that this procedure 
gives the accepted value for the energy per particle in the case of a pure Coulomb potential 
and a = 0. 

From the computations we found that the Laughlin state is favored for a pure Coulomb 
potential V(r) = e 2 /er. The actual potential between electrons in a real two-dimensional 
layer of finite thickness is softer, owing to the averaging of the charge density over the third 
dimension. This has been discussed in detail by numerous authors and the chice of potential 
depends on the shape of the well. We shall take the simple form of Zhang and Das Sarma. 
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They showed that V(r) = e 2 /e(r 2 + A 2 ) 1 / 2 is a reasonable approximation for a square 
quantum well and that A ~ 0.2t, where t is the layer thickness of the well. To understand 
the effect that this modification of the potential will have, note that the total energy is given 
by 

Ua = ^i (r 2 + A 2)l/ 2 ^- 1 ]- ( 3 ) 

We now plot the angle averaged correlation functions for different values of a in Fig. |2|. 
There is incipient solid order as a increases, in that there is much more tendency towards 
being able to identify shells of neighbors. At short distances, the correlation function is 
proportional to r 2 , just as that for the Wigner crystal is, and in contrast to the r 6 behavior 
of the Laughlin state. Overall, the the softer potential favors finite a: the correlation hole 
has a larger effective radius, even though it is not as "deep". The difference between the 
energy per particle for a BRS state with a = 1 and the Laughlin state for different values 
of A is plotted on Fig. |3|. There is thus a critical value A c at which the system undergoes a 
transition to finite a. We compute this to be A c = 4.1±1..5, which corresponds to a thickness 
of t = 1600A, when B = 10T. This transition is second order, unlike the transition to the 
Wigner solid, which is probably first order. It has recently been pointed out that changing 



t can induce the liquid-solid transition by a mechanism similar to that proposed here. flO 
The critical value of t for the m=l/3 density is similar to that computed here, suggesting 
that the energy balance between Laughlin, BRS, and crystalline states is a subtle one. We 
expect that the BRS state occupies a fairly narrow range of parameter space between the 
liquid and the crystal states, by analogy with hexatic phases. It is clear, in any case, that 
this range of thickness values is experimentally accessible. [|TTJ The energy balance between 
the BRS and Wigner crystal state is currently under investigation. |I2|] 

The correlations are oscillatory even to infinite distances in a true crystalline state. 
Inspection of Fig. |2| shows that the correlation function for the BRS state is flat at large dis- 
tances, demonstrating that long-range translational symmetry is not broken, and justifying 
the identication of these states as liquid crystal states. 



We may construct similar wavefunctions for m=5 and m=7. 

N 

^a(zi) = J[[(zi - Zj)(zi - Zj - a£)(zi - Zj + a£)(zi - Zj - (5£){zi - Zj + (5£)] exp(- ^ \zi\ 2 /At 2 ) 

i<j i 

(4) 

is an appropriate wavefunction for m=5 and 

N 

^a(zi) = ni(zi ~ z j){zi - Zj - a£)(zi - Zj + a£)(zi - Zj - (3£){zi - Zj + (3£) (5) 

i<j 

x ( Zi - Zj - 7 e){zi - Zj + 7 £)] exp(- 2 Iztf/Af) (6) 

i 

is an appropriate wavefunction for m=7. We have not yet investigated these wavefunctions 
for all values of (5 and 7. Particularly interesting cases are (5 — ia for m=5 and (5 = ua, 7 = 
u 2 a for m=7 if 00 is chosen as exp(27n/6). The polynomial parts of the wavefunctions for 
this parameter choice may also be written as fl- 2 ^ 4 — °- A ) (f° r m=5) and n^l- 26 — a& ) 
(for m=7), where The correlation functions for the m=5 wavefunction have a 

four-fold rotation axis and a six-fold rotation axis for the m=7 wavefunction. The incipient 
solid ordering again will stabilize these states as the well increases in thickness. For m=5 
and m=7 we calculate critical values A c (m = 5) = 2.9 ± 0.3 and A c (m = 7) = 2.1 ± 1.7, 
respectively. The latter is a state which resembles the hexatic state of two-dimensional 
fluids. The hexatic state, however, does not pick out particular directions in space while the 
m=7 state does. The m=5 state is somewhat similar to a biaxial nematic. 

The director n = (Re a, Im a) is the order parameter whose appearance signals the 
appearance of BRS. The states characterized by n and -n are identical. There is no indepen- 
dent inversion symmetry operation in our two-dimensional system and thus the transition to 
the low-symmetry phase is second-order, unlike the situation for ordinary three-dimensional 
nematic systems. Our Monte Carlo calculations of energy as a function of a confirm this 
picture. The Ginzburg-Landau energy is therefore 

F = A(T, t)n 2 + Bn 4 + K X (V ■ nf + K 2 (V x n) 2 , (7) 

where t is the thickness and T is the temperature. We expect a second-order transition 
when T=0 and our calculations have been carried out only at zero temperature. At finite 
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temperature thermal fluctuations are important. Depending on the experimental situation 
they may convert the transition to one of the Kosterlitz-Thouless type. 

The free emergy expression shows that twists of the order parameter are possible and 
lead to textures but also to low energy excited states. Nevertheless, these excited states 
do not involve density changes and the state as a whole is incompressible. We conclude 
that the FQHE still occurs in this gapless state. The BRS state does not appear to be 
a candidate for the Hall insulator phase. The quasiparticle and quasihole excitations are 
still gapped and their charges are the usual fractional ones. Their charge density profiles 
will have elliptical distortion. Since the projected oscillator strength f(k) and the projected 

— * — * 

static structure factor S(k) depend on the direction of k, the magnetoroton excitations with 

— * — * — * 

energy E(k) = f(k)/S(k) have dispersion which depends on the direction. 

From the experimental point of view, it appears that the chief difficulty in identifying the 
BRS states lies in distinguishing them from the Laughlin state. Correlation functions are 
anisotropic, but scattering experiments to test this are difficult to perform in two-dimensional 
systems. Tensor quantities such as the conductivity have a characteristic anisotropy (bire- 
fringence): <j xx {uj) 7^ a y y(uj), except at zero frequency, when a xx {u) = 0) = cr yy {uj = 0) = 0, 
as usual. Propagation of surface acoustic waves or measurements of the microwave conduc- 
tivity, perhaps with the simultaneous application of a current to eliminate domain effects, 
may be tools which can probe such an anisotropy, and test for the existence of the BRS 
states. 

We wish to thank A.V. Chubukov and M.B. Webb for useful discussions. This work was 
supported by the National Science Foundation under Grant No. DMR-9214739. 
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FIGURES 

FIG. 1. Typical configuration of the particles in the Monte Carlo simulation of the state given 
by Eqn. (1) with a = 3.2 The anisotropy of the correlations is clearly evident. 

FIG. 2. Angle averaged pair correlation function for the m=l/3 state for a = and a = 3.2. 
The a = 3.2 state shows incipient crystalline behavior at short distances and liquid-like behavior 
at long distances. 

FIG. 3. The difference between the energies per particle in a BRS state with a = 1 and the 
Laughlin state as a function of A, which is a measure of the well thickness. 
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